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Abstract. In this short note, we revisit the work of T. Tao and V. Vu on large non- 
hermitian random matrices with independent and identically distributed (i.i.d.) entries 
with mean zero and unit variance. We prove under weaker assumptions that the limit 
spectral distribution of sum and product of non-hermitian random matrices is universal. 
As a byproduct, we show that the generalized eigenvalues distribution of two independent 
matrices converges almost surely to the uniform measure on the Riemann sphere. 

1. Introduction 

We start with some usual definitions. We endow the space of probability measures on 
C with the topology of weak convergence: a sequence of probability measures (/i n )n>i 
converges weakly to /i is for any bounded continuous function / : C — > R, 



fdfin - J fd/j, 

converges to as n goes to infinity. In this note, we shall denote this convergence by 
fi n [j,. Similarly, for two sequences of probability measures (/%) n >i> (jJ>' n ) n >ij we wm 

n— >oo 

use /j, n — n' n 0, or say that \i n — fx' n tends weakly to 0, if 



J fdfin- J fdfi' n 

converges to for any bounded continuous function /. We will say that a measurable 
function / : C — > R is uniformly bounded for (fi n ) n ^i if 

limsup / \f\dfi n < oo. 

n— >oo J 

Finally, recall that a function / is uniformly integrable for (fJ> n )n^l if 

lim limsup / \f\dfi n = 0. 



The above definitions will also be used for probability measures on R + = [0, oo) and 
functions / : R+ -> R. 

The eigenvalues of an n x n complex matrix M are the roots in C of its characteristic 
polynomial. We label them Ai(M), . . . , A n (M) so that |Ai(M)| ^ ^ |A„(M)| ^ 0. We 
also denote by s\(M) ^ ••• ^ s n (M) ^ the singular values of M, defined for every 
1 < < n by Sfc(M) := \k{V MM*) where M* is the conjugate transpose of M. We 
define the empirical spectral measure and the empirical singular values measure as 

1 n 1 n 

^ M = ~Y1 6 ^ (M) and v M = - S Sk {M) . 

fc=l k=l 
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Note that [xm is a probability measure on C while vm is a probability measure on K + . 
The generalized eigenvalues of (M,N), two n x n complex matrices, are the zeros of the 
polynomial det(M — zN). If N is invertible, it is simply the eigenvalues of iV _1 M. 

Let (Xij)ij^i and (Yij)i j>i be independent i.i.d. complex random variables with mean 
and variance 1. Similarly, let (GV/)jj;si and (-Hy)ij>i be independent complex centered 
gaussian variables with variance 1, independent of (Xij, Yij). We consider the random 
matrices X n = (Xij)i^ij^ n , Y n = (Xij)i<,i,j^ni G n = (Gij)i<^ij^ n and H n = (Hij)i^ij^ n . 
For ease of notation, we will sometimes drop the subscript n. It is known that almost 
surely (a.s.) for n large enough, X is invertible (see the forthcoming Theorem II ip and 
then Hx~ 1 Y is a we U defined random probability measure on C. 

Now, let /i be the probability measure whose density with respect to the Lebesgue 
measure on C ~ I 2 is 

1 

7T(l + \z\ 2 ) 2 ' 

Through stereographic projection, [i is easily seen to be the uniform measure on the 
Riemann sphere. Haagerup and several authors afterwards have independently observed 
the following beautiful identity (see Krishnapur [15J, Rogers [19J and Forrester and Mays 

Lemma 1 (Spherical ensemble). For each integer n ^ 1, 

^G-^H = M- 

By reorganizing the results of Tao and Vu |21|, I22j. we will prove a universality result. 
Theorem 2 (Universality of generalized eigenvalues). Almost surely, 

n— >oo 

Applying once Lemma [T] and twice Theorem [21 we get 
Corollary 3 (Spherical law). Almost surely 

fi x -i Y ~» fJ.. 

n— >oo 

This statement was recently conjectured in |19[ IB]. More generally, our argument also 
leads to the following universality result for sums and products of random matrices. 

Theorem 4 (Universality of sum and product of random matrices). For every integer n, 
let M n , K n , L n be n x n complex matrices such that, for some a > 0, 

(i) x i — y x~ a is uniformly bounded for (vK n )n^i, and (vL„)n^i andx \— > x a is uniformly 
bounded for {v Mn )n>l, 

(ii) for almost all (a. a.) z G C, v K -i M l-i_k~ 1 l~ 1 z converges weakly to a probability 
measure v z . 

Then, almost surely, 

VM+KXL/VTi ^ ^ 

where [i depends only {y z )z^c- 

For M = K = L = I , the identity matrix, this statement gives the famous circular law 
theorem, that was established through a long sequence of partial results [171 13 CEH El El 
H EH El 13 [HI [211 122] • In this note, the steps of proof are elementary and they borrow all 
difficult technical statements from previously known results. Nevertheless, this theorem 
generalizes Theorem 1.18 in [22J in two directions. First, we have removed the assumption 
of uniformly bounded second moment for ^M+KXL/y/ni U K- 1 ML- 1 an d Vj^-i^-i. Secondly, 
it proves the convergence of the spectral measure. The explicit form of \i in terms of v z is 
quite complicated. It is given by the forthcoming equations ([21E]). This expression is not 
very easy to handle. However, following ideas developed in [19J or using tools from free 
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probability as in [20^ I12j . it should be possible to find more elegant formulas. For nice 
examples of limit spectral distributions, see e.g. [19]. It is interesting to notice that we 
may deal with non-centered variables (Xij) by including the average matrix of X/y/n into 
M, and recover [?]. Finally, as in [?], it is also possible by induction to apply Theorem [J] 
to product of independent copies of X (with the use of the forthcoming Theorem [5]) . 

2. Proof of Theorem [2] 

2.1. Replacement Principle. The following is an extension of Theorem 2.1 in [22]. The 
idea goes back essentially to Girko. 

Lemma 5 (Replacement principle). Let A n ,B n be nx n complex random matrices. Sup- 
pose that for a. a. z £ C, a.s. 

(i) va u -z - VB n -z tends weakly to 0, 

(ii) ln(-) is uniformly integrable for {y Ari -z)n>\ and (y Bn - z )n>\- 

Then a.s. p,A n -/^e„ tends weakly to 0. Moreover the same holds if we replace (i) by 

(i') fln(x)dvA n -z - f\n(x)dv Bn - z tends to 0. 

Proof. It is a straightforward adaptation of [3l Lemma A. 2]. □ 

Corollary 6. Let A n , B n , M n be n x n complex random matrices. Suppose that a.s. M n 
is invertible and for a. a. z £ C, a.s. 

(i) V A n -zM^ ~ V B n -zM^ tends Weakl V t0 °> 

(ii) ln(-) is uniformly integrable for (v An _ zM -i) n ^i and iy B zM - i) n ^i- 
Then a.s. HM n A n ~ ^M n B n tends weakly to 0. 
Proof. If M n is invertible, note that 

/ \n{x)dv M A z = — In | det(M n ^l n - z)\ = [ ln(x)du A _ zM -i + - In | det M n \. 
J n J n 

We may thus apply Lemma[^i')-(ii). Indeed, in the expression f \n(x)dvA n -z— J ^{x)dvB n -z-, 
the term ^ In | det M n \ cancels. □ 

2.2. Convergence of singular values. The following result is due to Dozier and Silver- 
stein. 

Theorem 7 (Convergence of singular values, [1]). Let (M n ) n ^i be a sequence of n x n 
complex matrices such that vm„ converges weakly to a probability measure v . Then a.s. 
v X n /Vn+M n converges weakly to a probability measure p which depends only on v. 

The measure p has an explicit characterization in terms of v. Its exact form is not 
relevant here. 

2.3. Uniform integrability. In order to use the replacement principle, it is necessary to 
prove the uniform integrability of ln(-) for some empirical singular values measures. This 
is achieved by proving that, for some /3 > 0, x i— > x~@ + x@ is uniformly bounded. 

Theorem 8 (Uniform integrability). Let (M n ) n ^i be a sequence ofnxn complex matrices, 
and assume that n->i" is uniformly bounded for (^M n )n>i for some a > 0. Then there 
exists (3 > such that a.s. x i— > x~@ + x 13 is uniformly bounded for (i / x„/ v / n+M„)n>l- 

In the remainder of the paper, the notation n>l means large enough n. We start with 
an elementary lemma. 

Lemma 9 (Large singular values). Almost surely, for n S> 1, 
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Proof. We have ^E7=l s i( X /V^) = ^ X *X = 4, £ 1<iJ<n \ X H | 2 , and the latter con- 
verges a.s. to 1 by the law of large number. □ 

Corollary 10. Let < a ^ 2 and Ze£ (M n ) n ^i 6e a sequence of n x n complex matrices 
such that x ^ x a is uniformly bounded for (fM„)n^l- Then, a.s. x h-» x a is uniformly 
bounded for (v Xn /y/n+M n )n>i- 

Proof. If M, iV are nx n complex matrices, from [1'6\ Theorem 3.3.16], for all 1 ^ i, j ^ n 
with l^i + j^n + 1, 

s i+j -i(M + TV) < Si(M) + Sj (N). 

Hence, 

s 2i (M + N) < s 2i -i(M + iV) < Si (M) + s;(iV). 
We deduce that for any non-decreasing function, / : M + — > M + and t > 0, 

/(x)aV A/+7V ^ 2 /" f(2x)dv M + 2 j f{2x)dv N) 

where we have used the inequality 

f(x + y)^f(2x) + f(2y). 

Now, in view of Lemma El we may apply the above inequality to f(x) = x a and deduce 
the statement. □ 

The above corollary settles the problem of the uniform integrability of ln(-) at +oo for 
v x/s/n+M- The uniform integrability at 0+ is a much more delicate matter. The next 
theorem is a deep result of Tao and Vu. 

Theorem 11 (Small singular values, [21, 22J). Let (M n ) n ^i be a sequence ofnxn complex 
invertible matrices such that x \— > x a is uniformly bounded for (^M„)n>l f or some a > 0. 
There exist c±,co > such that a.s. for n 3> 1, 

s n (X n /^ + M n ) ^ n- Cl . 

Moreover for i ^ n 1-7 with 7 = 0.01, a.s. for n>l, 

/\/n + M n ) ^ Co — . 

n 

Proof. The first statement is Theorem 2.1 in |21j and the second is contained in [22J (see 
the proof of Theorem 1.20 and observe that the statement of Proposition 5.1 remains 
unchanged if we consider a row of the matrix X n + ^/nM n ). □ 

Proof of Theorem^ By Corollary 1101 it is sufficient to prove that x 1— > x~@ is uniformly 
bounded for iyx/y/n+M) ano - som e j3 > 0. We have 

1 ™ 

limsup - Y) Sj'^/v'n + M) < 00, 
n n, 

1=1 

By Theorem llll we may a.s. write for n S> 1, 



i=l i=l i=Ln 1 -">'j4-l 



i=\ 

This last expression is uniformly bounded if < j3 < min(7/ci, 1). □ 
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2.4. End of proof of Theorem [2j If p is a probability measure on C\{0}, we define p as 
the pull-back measure of p under <fi : z \-¥ 1/z, for any Borel E in C\{0}, p(E) = p{(p^ 1 {E)). 
Obviously, if (p n )n>i is a sequence of probability measures on C\{0}, then p n converges 
weakly to p is equivalent to p n converges weakly to p. 

Note that by Theorem[51 a.s. for n>l, X n is invertible and x t— > x~@ + x@ is uniformly 
bounded for iyx n IJn)rC^\- Also, from the quarter circular law theorem, i^x n /y^i converges 
a.s. to a probability distribution with density 

IT 

(see Marchenko-Pastur theorem [161 E3 El])- From the independence of (Xij), (Yij), (Gij), (Hij), 
we may apply Corollary [6l Theorem [7J and Theorem [8] conditioned on (X^ ) to M n = 
zX n /^/n. We get a.s. 

Mx-iy - f J -x- 1 H ~^ 0. 

n— >oo 

By Theorem dH a.s. for n > 1, AT" 1 !? and G~ l H are invertible, it follows that 
Px- 1 H ~ f J 'G~ 1 H is equivalent to Px~ 1 H ~ f t G~ 1 H ~* 0. 

n— >oo n— >oo 

However since /xmjv = MiVM an d Am = /"a/- 1 ) we get 

Px~ x H ~ ^g~ 1 h ~* is equivalent to 1^h~ 1 x ~ P j H~ 1 g ~* 0- 

n— >oo n— >oo 

The right hand side holds by applying again, Corollary [6l Theorem [7] and Theorem [8l 

3. Proof of Theorem 3] 
3.1. Bounds on singular values. 

Lemma 12 (Singular values of sum and product). If M,N are n x n complex matrices, 
for any a > 0, 



J x a dv N ^J , 



X a dv M +N ^ 2 1+a ( J X a dv M + 

\ 1/2 / p \ 1/2 

X a dv M N ^21/ X 2a d^M / S 2a dl/JV 



Proof. The first statement was already treated in the proof of Corollary [TOj Also, from 
|13} Theorem 3.3.16], for all 1 ^ i, j ^ n with 1 ^ i + j ^ n + 1, 

Si+i _i(MAT) < s f (Af 

Hence, 

S2i (MiV) < s 2i _i(MAT) < Si (M) Si (iV). 

We deduce 



/ 



x a du MN <-Vs?(M)*?(JV). 

T7 ^ ^ 



n 
i=i 



We conclude by applying the Cauchy-Schwarz inequality. □ 

3.2. Logarithmic potential and Girko's hermitization method. We denote by T>'(C) 
the set of Schwartz distributions endowed with its usual convergence with respect to all 
infinitely differentiable functions with bounded support. Let 'P(C) be the set of probabil- 
ity measures on C which integrate In |-| in a neighborhood of infinity. For every p 6 V(C), 
the logarithmic potential of p on C is the function : C — > [— oo,+oo) defined for 
every z £ C by 

UJz) = / In \z - z'\ p(dz'), 
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(in classical potential theory, the definition is opposite in sign). Since In | - 1 is Lebesgue 
locally integrable on C, one can check by using the Fubini theorem that Uu is Lebesgue 
locally integrable on C. In particular, < oo a.e. (Lebesgue almost everywhere) and 
Ufj, G T>'(C). Since In |-| is the fundamental solution of the Laplace equation in C, we have, 
in V'(C), 

AU^ = TTfi, (1) 

where A is the Laplace differential operator on C is given by A = \{d^ + dy). 

We now state an alternative statement of Lemma [5] which is closer to Girko's original 
method, for a proof see O Lemma A. 2]. 

Lemma 13 (Girko's hermitization method). Let A ri be a n x n complex random matrix. 
Suppose that for a. a. z £ C, a.s. 

(i) VA n -z tends weakly to a probability measure v z on M+ ; 

(ii) ln(-) is uniformly integrable for [yA n -z)n^\- 

Then there exists a probability measure /x G V{C) such that a.s. 

(j) \iA n converges weakly to fj, 
(jj) for a.a. z G C, 



U^(z) = J ln(x) dv z . 



Moreover the same holds if we replace (i) by 
(i') J ln(x)duA n - z tends to f\n{x)dv z . 

Corollary 14. Let A n , K n , M n benx n complex random matrices. Suppose that a.s. K n 
is invertible and ln(-) is uniformly bounded for {vK n )n^i, and for a.a. z G C, a.s. 

(i) u a +k^ 1 (m -z) t en ds weakly to a probability measure v z , 

(ii) ln(-) is uniformly integrable for ( l/ A n +K' 1 (M n -z)) n ^ 1 - 

Then there exists a probability measure \x G V(C) such that a.s. 

(j) HM n +K n A n converges weakly to (j,, 
(jj) inV'{C), 

fi = —A / ln(x) dv z . 
^ J 

Proof. If K n is invertible, we write 

J \n(x)dv Mn +K n A n -z = j \n{x)dv An+K -i {Mn _ z) + ^ In | det K n \. 

By assumption, ^ln|deti^ n | = J ' \a{x)di>K n is a.s. bounded. We may thus consider any 
converging subsequence and apply Lemma [5ji')-(ii) together with (pQ). D 



3.3. End of proof of Theorem [4j We first notice that 

^M+KXL/^/n = ^LML-t+LKX/^/n- 

It is thus sufficient to prove that the right hand side converges. We set M = LML^ 1 and 
K = LK. Since K~ l {M - z) = K~ l ML~ l - K~ x L~ x z, we may apply Lemma [T2l and 
deduce that x h-» x a / A is uniformly bounded for (y^ _ z ))n>i- It only remains to invoke 

Theorem [5] and Theorem [7] applied to K~ 1 (M — z), and use Corollary [T5] for M + KX/ ' ^fn. 
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3.4. Explicit expression of the limit spectral measure. Let C + = {z G C : 3m(z) > 
0}, for a probability measure p on K, its Cauchy-Stieltjes transform is defined as, for all 

z G C+, 

1 



m p (z) = / dp(x). 

J x- z 

By Corollary [TJ] and Theorem in 2?'(C), 

li= ^-A J\n(x)dp z (x), (2) 

where for z G C, p z is a probability distribution on R + . From [I], for a. a. z £ C, p z has a 
Cauchy-Stieltjes transform that satisfies the integral equation: for all w G C+, 

. . f 2x(l + m (w)) , . . 
™*«0 = / 2 7TT 2 3 

J x l — (1 + m Pz (w)) z w 

where v z is as in Theorem UJ 



References 

1. Z. D. Bai, Circular law, Ann. Probab. 25 (1997), no. 1, 494-529. [2] 

2. Z. D. Bai and J. W. Silverstein, Spectral Analysis of Large Dimensional Random Matrices, Mathematics 
Monograph Series 2, Science Press, Beijing, 2006. [2] 

3. Ch. Bordenave, P. Caputo, and D. Chafai', Circular Law Theorem for Random Markov Matrices, 
preprint arXiv:0808. 1502v2, 2010.13111] 

4. B. Dozier and J. Silverstein, On the empirical distribution of eigenvalues of large dimensional 
information-plus-noise-type matrices, J. Multivariate Anal. 98 (2007), no. 4, 678—694. [3] [T] 

5. A. Edelman, The probability that a random real Gaussian matrix has k real eigenvalues, related distri- 
butions, and the circular law, J. Multivariate Anal. 60 (1997), no. 2, 203-232.(2] 

6. P. Forrester and A. Mays, Pfaffian point process for the gaussian real generalised eigenvalue problem, 
preprint arXiv : 0910 . 2531, 2009.0 

7. V. L. Girko, The circular law, Teor. Veroyatnost. i Primenen. 29 (1984), no. 4, 669-679.(2] 

8. , Strong circular law, Random Oper. Stochastic Equations 5 (1997), no. 2, 173-196.(2] 

9. , The circular law. Twenty years later. Ill, Random Oper. Stochastic Equations 13 (2005), 

no. 1, 53-109.(2] 

10. I. Y. Goldsheid and B. A. Khoruzhenko, The Thouless formula for random non-Hermitian Jacobi matri- 
ces, Israel J. Math. 148 (2005), 331-346, Probability in mathematics. MR MR2191234 (2006k:47082) 

m 

11. F. Gotze and A. Tikhomirov, The Circular Law for Random Matrices, preprint to appear in the Annals 
of Probability arXiv : math/07093995 [math . PR] , 2010. d 

12. A. Guionnet, M. Krishnapur, and O. Zeitouni, The single ring theorem, preprint arXiv : 0909 . 2214, 
2009. 

13. R. Horn and C. Johnson, Topics in matrix analysis, Cambridge University Press, Cambridge, 1991.2] 

m 

14. C.-R. Hwang, A brief survey on the spectral radius and the spectral distribution of large random matrices 
with i.i.d. entries, Random matrices and their applications (Brunswick, Maine, 1984), Contemp. Math., 
vol. 50, Amer. Math. Soc, Providence, RI, 1986, pp. 145-152.(2] 

15. M. Krishnapur, From random matrices to random analytic functions, Ann. Probab. 37 (2009), no. 1, 
314-346. (3 

16. V.A. Marchenko and L.A. Pastur, The distribution of eigenvalues in sertain sets of random matrices, 
Mat. Sb. 72 (1967), 507-536.(5] 

17. M. L. Mehta, Random matrices and the statistical theory of energy levels, Academic Press, New York, 
1967. [2] 

18. G.M. Pan and W. Zhou, Circular law, extreme singular values and potential theory, J. Multivar. Anal. 
101 (2010), no. 3, 645-656.(2] 

19. T. Rogers, Universal sum and product rules for random matrices, J. Math. Phys. 51 (2010), no. 093304. 

(3 El E 

20. P. Sniady, Random regularization of Brown spectral measure, J. Funct. Anal. 193 (2002), no. 2, 291- 
313.0 

21. T. Tao and V. Vu, Random matrices: the circular law, Commun. Contemp. Math. 10 (2008), no. 2, 
261-307.(211] 



s 



CHARLES BORDENAVE 



22. , Random matrices: Universality of ESDs and the circular law, preprint to appear in the Annals 

of Probability arXiv : 0807 . 4898 [math . PR] , 2010. H El HI 

23. K. W. Wachter, The strong limits of random matrix spectra for sample matrices of independent ele- 
ments, Ann. Probability 6 (1978), no. 1, 1-18. [5] 

24. Y. Q. Yin, Limiting spectral distribution for a class of random matrices, J. Multivariate Anal. 20 
(1986), no. 1, 50-68. M 

Acknowledgment 

The author is indebted to Tim Rogers for pointing reference [19] which has initiated this 
work, and thanks Djalil Chafai and Manjunath Krishnapur for sharing their enthusiasm 
on non-hermitian random matrices. 

(Ch. Bordenave) IMT UMR 5219 CNRS and Universite Paul-Sabatier Toulouse III, France 
E-mail address: charles .bordenave (at)math.univ-toulouse . fr 
URL: http: //www .math. univ-toulouse . f r/ "bordenave/ 



